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Abstract. We introduce a multi-parameter generalization of the Lambda-determinant of 
Robbins and Rumsey, based on the cluster algebra with coefficients attached to a T-system 
recurrence. We express the result as a weighted sum over alternating sign matrices. 
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I. Introduction 

The so-called Lambda-determinant was introduced by Robbins and Rumsey [l4j as a 
natural generalization of the ordinary determinant, via a one-parameter deformation of the 
Dodgson condensation algorithm [7] that expresses the determinant of any n x n matrix 
in terms of minors of sizes n — 1 and n — 2. This produces, for each n x n matrix A, a 
Laurent polynomial of its entries, involving only monomials with powers of ±1, coded by 
Alternating Sign Matrices (ASM) of same size n. Robbins and Rumsey were able to write 
a compact formula for the Lambda-determinant of any matrix A, as a sum over the ASMs 
of same size with explicit coefficients [T5] (see also [2] for a lively account of the discovery 
of ASMs). ASMs are known to be in bijection with configurations of the Six Vertex (6V) 
model on a square nxn grid, with special Domain Wall Boundary Conditions (DWBC) |10j . 
The latter are obtained by choosing an orientation of the edges of the underlying square 
lattice, such that each vertex of the grid has two incoming and two outgoing adjacent edges. 
The resulting six local configurations read as follows: 



i | | 1 1 v 

6V ^^^^r^^- 

a 1 a 2 b 1 b 2 c 1 c 2 


6V 4^ 
DWBC 3 : 


• 
» 
» 


* T*" 

• m <t 


ASM 00001-1 









The boundary condition (DWBC) imposes that horizontal external edges point towards 
the grid, while vertical external edges point out of the grid. The bijection between the 
6V configurations and ASMs goes as follows: the a and b type vertices correspond to the 
entries 0, while the c\ type vertex corresponds to the entry 1, and the C2 type vertex to the 
entry —1. Conversely, for any nxn ASM B, there is a unique 6V-DWBC configuration on 
the nxn square grid. For any fixed ASM B, we denote by C(B)ij the configuration of the 
vertex in the corresponding 6V-DWBC model, with Q(B)ij G {aj., a2, b±, 62, ci, C2}. 

In this note, we define an inhomogeneous, multi-parameter generalization of the Lambda- 
determinant of an n x n matrix A. This is done by deforming the Dodgson condensation 
algorithm [7] by use of two sets of parameters A a ,/z a , a G Z. More precisely: 

Definition 1.1. Let A = (ajj)jj e / ; i = {1, 2, n}, and sequences A = (\ a ) a &, H = 
(^ a )a& of fixed parameters. The generalized Lambda- determinant of A, denoted \A\\ )fJ , is 
defined inductively by the following modified Dodgson condensation algorithm. Let Af^'"' 3 * 
denote the submatrix of A obtained by erasing rows i±,...,ik and columns and s 

be the shift operator acting on sequences as (s\) a = \ a +i and (s _1 A) a = A a _i and similarly 
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on fi. We have: 

I . I A^O 1^4-1 |sA,s/x l^nls-iA,*- 1 /^ + A() | ^4™ | s" 1 A,s/i | | sX^' 1 /! 

T/ws recursion relation on the size of the matrix, together with the initial condition that a 
0x0 matrix has generalized Lambda- determinant 1 and a 1 x 1 matrix (a) has generalized 
Lambda- determinant a determines \A\\^ completely. 

This generalization of the Lambda-determinant is connected to the expressions found 
in [15] for solutions of the octahedron recurrence as partition functions of domino tilings 
of Aztec diamonds. The main result of the present paper is the following closed formula, 
which holds whenever the definition makes sense: 

Theorem 1.2. For fixed parameters A a ,/i a; a G {— n + 2, — n + 3, ...,n — 3,n — 2} the 

generalized Lambda-determinant of a matrix A = (djj)i.je/, / = {1, 2, n}, is expressed 
as follows: 

n 

(1-2) \A\^= Yl n^CASjA^) 

ASM B i,j=l 

where the sum extends over all the n x n Alternating Sign Matrices B = (bij) it j e j, and the 
weights Wij(A, B, A; yu) are defined as: 

{Xj^i if G(B) id = ai 

/^j+j-n-i if ^\B)i j = b\ 

(V.+W.-i) ife(B)y = c 2 - 

1 otherwise 

The paper is organized as follows. In Section 2, we recall Robbins and Rumsey's 
original definition of the Lambda-determinant, and express it as a solution of the T- 
system/octahedron recurrence with one coefficient. In Section 3, we define the generalized 
Lambda-determinant as a solution of a T-system with inhomogeneous coefficients. The 
latter is shown to be part of a cluster algebra of infinite rank, and as such to possess the 
Laurent property. Using the T-system relation, we present a non-homogeneous example 
leading to simple product formulas for the generalized Lambda-determinant of the matrix 
with entries a it j = 1 for all i, j. Section 4 is devoted to the actual computation of the gen- 
eralized Lambda-determinant, by solving the T-system with inhomogeneous coefficients. 
This is done via a matrix representation of the T-system relation that leads to a general- 
ization of the solution of [B]. In Section 5 the solution is rephrased in terms of networks, 
namely of paths on directed graphs. The latter are finally mapped onto configurations of 
the 6V-DWBC model, leading to the proof of Theorem 11.21 Section 6 is devoted to further 
properties of the generalized Lambda-determinant, an expression purely in terms of ASMs, 



4 



PHILIPPE DI FRANCESCO 



some explicit examples, and the complete solution of the T-system with inhomogeneous 
coefficients, and finally some concluding remarks. 

Acknowledgments. We thank R. Kedem for many useful discussions. We also thank 
the Mathematical Science Research Institute in Berkeley, CA and the organizers of the 
semester "Cluster Algebras" for hospitality while this work was completed. This work is 
partially supported by the CNRS PICS program. 

2. The classical Lambda-determinant 

2.1. Desnanot-Jacobi identity and determinant. Let I = {1, 2, n}. For any matrix 
M = (mij)ijizj we use the same notations as in Def Jl.ll For any given k + 1 x k + 1 matrix 
M, we have the celebrated identity: 

(2.1) |M| x \M\£\ = |M™| x \M\\ - \Mf \ x \Ml\ 

This gives rise to the Dodgson condensation algorithm for effective computing of determi- 
nants, as the identity may be used as a closed recursion relation on the size of the matrix, 
allowing for computing its determinant from the initial data of determinants of matrices of 
size and 1 (equal respectively to 1 and the single matrix element). 

More formally, we may recast the algorithm using the so-called T-system relation: 

(2.2) T*,j,fc+iTij,fc-i = Ti,j+i,fcT;,j-i,fc — Ts+i,j,fcTi-ij',fc 

for any i, j,k G Z with fixed parity of i + j + k. 

Let A = (a>i,j)i,jei be a fixed n x n matrix. Together with the initial data: 

Te,m,o — 1 (£,m £ "Z; £ + m = n mod 2) 

(2.3) Tj^i = d j-i+n+i i+j+n+i (z, j £ Z; i + j = n + 1 mod 2; \i\ + \j\ < n — 1) , 

the solution of the T-system (12.21) satisfies: 

(2.4) T 0>0>n = det(A) 

2.2. T-system with a coefficient and Lambda-determinant. Robbins and Rumsey 
have introduced the Lambda-determinant by applying the Dodgson algorithm to the follow- 
ing modified Desnanot-Jacobi identity. Let A G C* be a fixed parameter and A = (aij)ij^i 
be a fixed n x n matrix. 

Definition 2.1. The Lambda-determinant of A, denoted \A\\ is defined as the solution 
\A\\ = T 0i o,n of the following T-system with a coefficient 

(2.5) T i j^ +1 T i j t i i ,_i = Tij + i^Tij_i^ + \ T i+ ij^Ti_ijk 
and subject to the initial conditions (12. 3p . 

Note that for A = —1, we recover the usual determinant: \A\-i = det(A). 
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Example 2.2. The generalized Lambda- determinant of a 3 x 3 matrix A reads: 

bdfh 



d e f 
g h k 



a b c 

= A 3 ceg + X 2 cdh + X 2 bfg + aek + X bdk + X afh + X(X + 1) 



Example 2.3. The Lambda- determinant of the Vandermonde matrix A = (a{ reads 

\A\ X = \ { (aj + \di) 

l<i<j<n 

2.3. The Robbins-Rumsey formula and Alternating Sign Matrices. A remarkable 
property of the definition of the Lambda-determinant, is that it produces a Laurent poly- 
nomial of the matrix elements a^j, which is a sum over monomials with only powers ±1. 
These powers are coded bynxn so-called Alternating Sign Matrices (ASM), namely ma- 
trices B = (bij)i.j & i with entries b{j 6 {0,1,-1}, with non-negative partial row sums: 
Z/j=i > 0, k = 1,2, ...,n — 1, i 6 /; and with row sums equal to one: YTj=\^i,j = 1> 
i 6 I. More precisely, Robbins and Rumsey found the following formula for the Lambda- 
determinant: 

(2.6) \A\ X = A^ll + A-f^n^ 

ASMB 

where 1)b is the total number of entries in B that are equal to —1, and 

Inv(.B) = ^2 b k,pbi, m 



k<l 



is the generalized inversion number of B. 

Remark 2.4. Note that for X = —1 the sum truncates to only the contribution of ASMs 
with no —1 entry, which are the permutation matrices P, and Inv(P) is the usual inversion 
number of the permutation matrix P, so that (I2.6P reduces to the usual determinant formula. 

Remark 2.5. From eg. fl2.6p . we see that the Lambda- determinant is well defined for any 
matrix A with non-vanishing entries. More precisely, as the matrix elements b^j = —1 of 
ASMs may only occur for 2 <i,j < n — 1, the Lambda- determinant is well-defined for any 
n x n matrix A such that a^j ^ for all 2 < i, j < n — 1. 



3. AN INHOMOGENEOUS GENERALIZATION OF THE LAMBDA-DETERMINANT 

3.1. T-system and Cluster Algebra. The unrestricted A^ T-system (I2.2p is known 
to be part of an infinite rank cluster algebra with coefficients [9]. The T-system in this 
form (with Tij^ defined only for i + j + k = 1 mod 2), was shown to be a particular 
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mutation in an infinite rank cluster algebra with coefficients, with cluster variables of the 
form (Tij jk _i, Tii j\k)i ,j,i',fez an d non-vanishing fixed coefficients (A;)j G z It reads: 



(3.1) 



As a consequence of the cluster algebra Laurent property [9], the solution Tij^. of this 
equation may be expressed as a Laurent polynomial of any choice of admissible initial 
values, with coefficients in Z[(A;)j e z]- The "coefficients" Aj are inhomogeneous, as they 
now depend on the value of i. In particular, choosing homogeneous coefficients A« = A 
independent of i, we recover the Laurent property observed for the Lambda-determinant 
of Robbins and Rumsey. 

As it turns out, we may consider an even more general inhomogeneous coefficient T- 
system equation: 



(3.2) 



-^i j,fc+l-^i,j,fe— 1 fJij -^ijj+ljA-^ij— l,k Aj Ti+x^ kTi— l,j,k 



for any fixed non-vanishing coefficients (Xi)iez and We have the following 



Theorem 3.1. The equation (13. 2 p is a particular mutation in a cluster algebra A with 
coefficents. The initial seed is given by the cluster Xq = ((T^cOijez, {Ti',j',i)i',j'ez) with 
i + j = 1 mod 2 and i' + f = mod 2 as usual, and the coefficients ((A a ) ag ^, (nb)bez) , and 

( B\ 
-B 
Lq Li 
\M M x j 
where: 



the extended exchange matrix B with the following infinite block form: B 



{Lo)a;i,j — &i,a (-^l)a;i',j' 



{M )b;i,j 



-5, 



(Mi 



l)b;i',j' 



-- 6f,b 



Proof. We must show that the cluster X k = ((T iJtk ) itjeZ , (Ti'j',fc+i)i',i'ez) , where T i>j)k is 
the solution of (13. 2p with prescribed values of ((T^oKjez, (^V^i^/j/gz), is a cluster in the 
cluster algebra A. First let us represent the quiver Q coded by the extended exchange 
matrix Bq. We represent by a filled • (resp. empty o) circle the vertices (i, j) such that 
i + j = 1 mod 2 (resp. (i', j') such that i' + j' = mod 2), and by a cross <8> (resp. square 
□) the vertices (a) indexing A a (resp. (b) indexing /i^) of Qq. We have the following local 
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structure, around respectively vertices • and o: 




Note that the vertices (a = i) indexing Aj are connected to all j £ Z while vertices 

(b — j) indexing fij are connected to all i G Z, with edges of alternating orientations. 

We may reach the cluster X\ from the initial cluster X by a compound mutation, obtained 
by mutating all the vertices •. The resulting quiver Q\ is identical to Qo with all edge 
orientations flipped, namely Bi = —B , as illustrated below for mutations in a quadrant 
around a • vertex: 




(Recall that a mutation at a vertex v of a quiver flips the orientations of all incident edges, 
creates "shortcut" edges u — > w for any length 2 path u — > v — > w before mutation, and 
any 2- loop v — > w — > v thus created must be eliminated.) 

The cluster X 2 is then reached by mutating all vertices o, and has the exchange matrix 
B 2 = B , etc. □ 

As a consequence of the theorem, we have the Laurent property: the solution T^j^ of the 
system ( 13. 2ft for prescribed values of ((^ij,o)tjeZ) (^i'j',i)t',j'ez), is a Laurent polynomial of 
these values, with coefficients in Z[(\ a ) aeZ , (fj, b ) beZ \. 

3.2. Generalized Lambda-determinant: definition and properties. In view of the 
results of previous section, we may recast the Definition 11.11 into the following: 

Lemma 3.2. Given a square matrix A = (ajj)ij e j, and some parameters A a ,// Q; a G J = 
{ — (n — 2), — (n — 3), — 1, 0, 1, n — 3, n — 2}, the inhomogeneous Lambda- determinant of 
A is the solution To 5 o,n of the inhomogeneous T -system ( I3.2p . with initial conditions ( 12. 3p . 
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As a consequence of Theorem 13.11 we know that |^4|a ; ^ is a Laurent polynomial of the 
entries a^j of A, with coefficients in Z[(A a ) ag j, (fJL a ) a ej\- 



d e f 
g h k 



Example 3.3. The generalized Lambda- determinant of a 3 x 3 matrix A reads: 
a b c 

= AiAoA_i ceg + Ao/^oAi cdh + Ao/ioA_i bfg 

bdf h 

+yUi/i yU-i aek + X /^o^i bdk + X /j ^-i afh + X n (X + Ho) 

e 

Example 3.4. Let X a = A and fi a = fi, for all a, and A be the Vandermonde matrix 
A = (a{~ 1 )ij ( zj. Then we have: 

\A\\. )ft = | j (Xai + fxaj) 

l<i<j<n 

The generalized Lambda-determinant reduces to the original Robbins-Rumsey Lambda- 
determinant for the choice of parameters A a = A and = 1 for all a. 

The generalized Lambda-determinant satisfies a number of properties, inherited from the 
symmetries of the T-system. We have the following 

Proposition 3.5. Let a, r be the following transformations of the matrix A = (a it j)ij e j, 
corresponding respectively to quarter-turn rotation, and diagonal reflection: 

Cr (A)i,j = a n+l-j,i = a j,i 

Then with the following action on parameters: 

(3.3) fo-(X) =n ^ fr(A) =A_ Q 

[ cr(/x) Q = A„ [ r(fi) a = n a 

we have: \<p{A)\ v ^),<pOi) = I^|a, m for ip = a,r. 

Proof. To compute |^(^4) | ¥ ,(a) ) ^(/^) , we use the following T-system relation 

with Xf = <f(X)i and /xj = y?(/i)j, together with the initial data 

= 1 (i + j — n mod 2) T^ fl = (p(A) j-»+n+i i+i+n+i (i + j — n + 1 mod 2) 



rptp rplf _ ..Vrptp rplfi i \lfirpip rp(fi 

± i,j,k+l ± i,j,k-l ~ Hj ± i,j+l,k ± i,j-l,k T A i 1 



We wish to compare the solution Tq Qti with that, T 0>0)n , of the T-system (13. 2p subject to 
the initial conditions (I2.3p . For ip = cr, r, we have respectively: 



Tfj 1 — °~(A) j-i + n+l i+j + n + l — a -i-j + n + 1 -i+j+n + 1 — Tj _j 

2 ' 2 2 ' 2 

TJ; 1 = T (A) j-i+n+l i+j + n + 1 = a i+j + n+1 -i+j + n + 1 = T_j 1 

'J' 9 ' 9 9 ! 9 
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From the obvious symmetries of the T-system, and with the respective choice of parameters 
( 13.30 . we deduce that 

^i,j,k = Tj,-i,k and T[ jk = T_ ijjjk 

for all i,j, k such that i + j + k = n mod 2. The proposition follows for i = j = and 
k = n. □ 

3.3. A non-homogeneous example. In this section we study a non-trivial example of 
an inhomogeneous generalized Lambda-determinant that may be of interest as a statistical 
model. We choose the matrix A = (ai t j)ij e j with entries 

(3.4) a M = l (ijel) 

We pick coefficients with an explicit dependence on their index, namely 

(3.5) X a = q a fi b = q b (a, b G Z) 
for some fixed parameter q G C*. We have the following 

Theorem 3.6. The generalized Lambda- determinant of the matrix A (13. 4p and with the 
parameters A, fi of (13.51) reads: 

L"2"J n—2m 

i^=n n 

m=l j=2m—n 

To prove the theorem, let us actually solve the more general T-system (I3.2p with coeffi- 
cients A a ,/Xfe as in ( 13.50 . and with initial conditions ( 12.30 . We have: 

Lemma 3.7. The solution of the inhomogeneous T-system (13.21) with initial data Tj Ji0 = 
Tij,i = 1 for all i,j G Z reads: 

l^P^-i k-\i-j -2m 

T Mtk = n n ( i+ ^ a ) 

m=l a=2m-k+\i— j\ 
k— \i— j\+2m— 2 

(3.6) x n (l + g a ) 

m=l a=m 

Proof. By uniqueness of the solution for the given initial data, we simply have to check 
that the above satisfies both the initial condition Tj^o = Tj J; i = 1 and the T-system with 
coefficients. The condition T iji0 — = 1 is clear from the formula ( 13.60 . Next, we note 
that the expression (13.60 satisfies: 

q 7r +l ' kT r' l,k = ^ iiJt*,keZ + ) 

Q 1 i-\-l,j,k 1 i-l,j,k 
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Analogously, we compute: 



9 J-i+l,j,kJ-i-l,j,k 



>0j 



and therefore (13. 2p follows, the lemma is proved. □ 



Theorem 13.61 follows by taking % = j = 0, k = n in the expression (I3.6P above. Note that 
the result of Theorem 13.61 reduces in the homogeneous limit q — > 1 to: 



n(n— 1) 

2o,0,n = 2 2 



in agreement with the result of Example 13.41 for = 1 for all i and A = \x = 1 . 



4. Computing the generalized Lambda-determinant 

We first write a solution of the T-system with inhomogeneous coefficients, based on a 
matrix representation generalizing the solution of [6]. 



4.1. Inhomogeneous T-system solution I: V and U matrices, definition and prop- 
erties. Let us consider the square lattice with vertex set Z 2 , and its elementary triangu- 
lations, obtained by picking either of the two possible diagonal edges in each square face. 
We consider a pair of such triangles sharing an horizontal edge of the lattice as a gener- 
alized rhombus, and we restricted ourselves to bicolored triangulations such that exactly 
one of the two triangles in each generalized rhombus is grey, the other one being white. 
We moreover attach variables to the vertices of the lattice This gives rise to eight possible 
generalized rhombi. In analogy with the solution of the A r T-system [3], we associate the 
following 2x2 matrices to each of these generalized rhombi, with entries Laurent monomi- 
als of the variables at the three vertices adjacent to the grey triangle, and also depending 
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on fixed coefficients A,//: 



V(d, a, b; A, ji) = 




d d 



V'(d,a,b; A, ji) = 



U(a,b,c; A, fi) = 



U' (a, 6, c; A, fi) = 




Note that we have represented a thicker blue horizontal (resp. red vertical) edge to indicate 
the presence of the parameter \x (resp. A) in the corresponding matrices. Next, we introduce 
a graphical calculus, by associating to pictures some matrix products. We read pictures 
from left to right, and take the product of matrices in the same order. 

With this rule, the above matrices satisfy the following property, easily checked by direct 
calculation. 

Lemma 4.1. 




iff xx 1 = fiab + Xcd 
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We use the standard embeddings of any 2x2 matrix M into SL n denoted Mj, % 
1, 2, n — 1 with entries 



(A*<) r , a 



5 r>s if r g" {i, z + 1} or s g" {i, i + 1} 



M, 



r, s 



otherwise 



In this embedding, it is clear that for any 2x2 matrices M, P Mi and Pj commute unless 
j — i G {0, ±1}. We have also the following exchange properties, easily checked directly: 

Lemma 4.2. The matrices U{,Vj satisfy: 

Ui(a, b, c; A, ii)V i+ i(b, c, d; A, /i) = V- +1 (a,c,d; A, fi)U-(a,b,d; A, /x) 





where again the pictures are read from left to right and the direction of the diagonal of a 
square indicates which generalized rhombus is to the left of the other, and we have repre- 
sented the 4 up/down pointing white triangles in-between the two possible configurations of 
half-squares each of them may be in (total of 16 configurations here). We also have the 
commutation relations: 

[Vi^a^.U^eJ)] = [V:(c,a,b),U,(d,e,f)} = [^(c, a, b), U'^d, e, /)] 

= [Vj(c,a,b),U' j (d,e,f)] = 

The relations of Lemmas 14.11 and 14.21 are all we need to construct the solutions of the 
unrestricted T-system with coefficients, expressed in terms of products of U and V 
matrices. The relation of Lemma 14.11 allows to implement the T-system evolution on 
the variable in the center x — > x'. The relations of Lemma 14.21 allow for rearranging 
triangulations so as to be able to apply the former. In particular, when representing white 
and grey triangulations, we may omit the diagonal of each square of a uniform color, as 
both choices of diagonal lead to the same matrix product. So we will typically consider 
tilings of domains of Z 2 by means of grey and white unit squares (2 tiles), and grey and 
white triangles equal to half a unit square, cut by one diagonal (4 tiles) such that any edge 
is common to two tiles of different colors. We call such a tiling a square-triangle tiling. 

4.2. Inhomogeneous T-system solution II: a determinant formula. We consider 
the tilted square domain DJi,j) C Z 2 defined by 



Dk(j, i) = {{x, y) G Z 2 such that \x — j\ + \y — i\ < k — 1} 



,j G Z; k G Z^ 



AN INHOMOGENEOUS LAMBDA-DETERMINANT 13 

We also consider two square-triangle tilings of Dk(j,i), denoted by 9 m in{k) and 9 max (k) 
defined as follows. O m i n (k) is the intersection of a checkerboard of alternating grey and 
white squares with D k (j,i) such that the NW border of D k (j,i) is adjacent only to white 
triangles (likewise, the SE border is adjacent to only white triangles, while the SW and 
NE borders are adjacent only to black triangles). 9 max (k)is the only square-triangle tiling 
without squares having opposite color boundary assignments (NW,SE grey, SW,NE white). 
Here are examples of 9 m in(k) and 9 max (k) for k = 5: 



9 rain (5) 



9 max (5) 




To attach matrices to these tilings, let us also include thick (red or blue) edges like in 
the representations of U, V, U', V . We fix the ambiguity on grey squares of 9 m i n (k) (from 
Lemma [4.21) by picking the l.h.s. representation. This gives two decorated tilings 9 m i n {k) 
and 9 max (k) which for k = 5 read: 



9 rain (5) 




9 max (5) 




We finally attach variables to the vertices and coefficients to the thickened edges of the 
tilings as follows. On 9 m in{k) we assign variables t%j to the vertices G Dk(0,0), 

and coefficients Aj to thickened vertical (red) edges with top vertex of the form (£,i), fij 
to thickened horizontal (blue) edges with right vertex of the form (j,m), £, m G Z. On 
9 m ax{k) we assign variables u^j to the vertices G -D fe (0,0), and coefficients Aj to 

thickened vertical (red) edges with top vertex of the form (£,i), fij to thickened horizontal 
(blue) edges with right vertex of the form (j, m), £, m G Z. For completeness, we assign the 
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value /i = 1 (resp. A = 1) to coefficients corresponding to non-thickened horizontal (resp. 
vertical) edges. 

We finally attach matrices Q m i n (k) and Q ma x{k) to the two tilings. These are defined 
as the products of Z7j, Vi, U-, V( matrices embedded into SL 2n -2, coded respectively by the 
two decorated tilings with assigned vertex and thickened edge values. 

Example 4.3. For k = 3, we have: 

(3) = V 2 (t 

x ^2(^0,-1 ? £o,o> £1,-1, 1) 1) ^3(^0,0) £1,-1, £i,o! Ai, /io) 

X[/i(t_L,Q, £-1,1, £o,05 1) l)^2(£-l,l) £o,Q) £o,i; Ao, /il) 

x ^3(£i,0) £i,i) £2,0; 1) l)^2(£o,l) £o,2) £1,1; 1, 1) 
corresponding to the following variables and coefficients: 



i 




-2-10 12 



Likewise, we have: 

@maz(3) = C/2(wo,-2, «Q,-1, «l,-i; l,/i-l)^l(W-l,-l, W-1,0, W0,-i; 1) 1) 

-l) «i,o, ^2,0; 1) ^0)^2(^0,-1, «o,o, «i,o; 1, yUo) 
xV2(«_i,o, wo,o, «o,i; A , l)V((u^ 2 ,o, W-1,0, W-1,1! A_i, 1) 

Xl^(wo,lWl,Q, «l,i; Al,//!)^'^-!,!,^,!, Wo,2| A , 1) 
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corresponding to the following variables and coefficients: 




M 1 1 1 1 - i 

-2-10 12 

We have: 

Theorem 4.4. Let n > 1, and let us pick variables ttjj as: 
(4-1) Uij = T ij>HiHi | (j>')eD n (0,0) 

where T^^, for i + j + k = n mod 2, is the solution of the T -system with coefficients ( 13. 21) 
subject to the initial conditions 

(4-2) T i>ji£ .. n = t itj (j,i)eD n (0,0) 

where £i t j n — (i + j + nmod2) 6 {0, 1}. Then we have the matrix identity 



Proof. The identity is proved by showing that we may transform the minimal tiling 6 m i n (n) 
with variables tij into the maximal tiling 8 max (n) with variables Uij by successive applica- 
tions of the Lemmas 14.11 and 14.21 This is illustrated below for the case n = 3. 




(a) (b) (c) (d) (e) 



Starting from 6 m i n (n) (a), let us first apply Lemma fl~2l to flip the diagonals of all the white 
squares (b). Then apply Lemma I4TT1 to update all the values at vertices (j, i) G -D n (0,0) 
such that i + j = n mod 2 (marked as empty circles on (b)). Due to the update relation 
xx' = fiab + Xcd, and comparing with the T-system evolution (13.21) . we see that the updated 
values are Tij,2 (filled vertices on (c)). We then use Lemma [4.21 to flip the diagonals of all 
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the remaining squares (d), then update the values at vertices G -D„_i(0,0) such that 
i + j = n — 1 mod 2 (empty circle on (d)), yielding updated values Tj Jj3 . The procedure 
is iterated until there are no more squares. In the final tiling configuration, the boundary 
triangles have opposite colors from those of 9 min (n), the tiling is therefore 6 max {n). Note 
that the last (central) updated value is m ,o — To,o,n- As both Lemmas 14.11 and 14.21 leave 
the result of the corresponding matrix products unchanged, the theorem follows. □ 

Corollary 4.5. Let a m i n {n), a max (n) denote respectively the n — lxn — 1 principal minors 
of Q m in(n) and Q max (n), then we have: 



[n = a m „An) 



Our last task is to relate the central value T 0i0i n = «o,o to the matrix 0™„(n). We have 
the following: 



Theorem 4.6. Notations are as in Theorem 4-4' The solution T 0) o,n of the T system with 
coefficients (13.21) is expressed in terms of its initial data as: 

-1 o 



^0,0,n — O m i n (n) J~| tji_ n _j ] [ t>i,n- 



-l+i 

i=2—n i=2—n 

Proof. Let us compute a max {n). We note that, as Ui, U[ are lower triangular and Vi, V( are 
upper triangular, the 2n — 2 x In — 2 matrix O ma x{n) is the product of a lower triangular L 
matrix by an upper one IX, corresponding respectively to the left and right halves of 6 max {n). 
More precisely, all matrices in the SW corner are of Ui type (i < n — 2), all matrices in the 
NW corner are of U[ type {n — 1 < i < 2n — 3), all matrices in the SE corner are of V( type 
(z < n — 1) and all matrices in the NE corner are of Vi type (n < i < 2n — 3). 

Let M max (n) be the truncation of ma z(^) to its n — 1 first rows and columns. From the 
triangularity property, we have M max (n) = £u where i and u are the truncations of £ and IX 
to their n — 1 first rows and columns, so that a max {n) = det(£) det(w). By the triangularity 
property, and the above remark on corners, det(£) is the product over all diagonal matrix 
elements of the Ui in Q m ax{n), while det(w) is the product over all diagonal matrix elements 
of the V( in Q m ax{n). These in turn have a very simple interpretation in terms of the tilings: 
the diagonal matrix elements of Ui (i < n — 2) are the ratios (left vertex value) / (right 
vertex value) along the horizontal edges of the tiling (here in the strict SW corner, i.e. 
with 2 — n < i < and 1 — n < j < 0. The product over these is telescopic and leaves us 
with only the ratios of (leftmost vertex value) /(rightmost vertex value) =Uj i j+i_ n /^j,o along 
the row i. Analogously, the diagonal matrix elements of V( (i < n — 1) have the same 
interpretation, and we get the telescopic products «i,o/wi,n-i+i along the row i. Collecting 
both products, we get 

-10 -10 

Omax (n) = det(£) det(ti) = FT JJ = «0,0 IT Ui,l-n-i IT Kn-l+i 

UiQ Ui n -i + i 
i=l— n ' i=l— n ' i=2— n i=l—n 
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Finally, we note that the SW and SE boundary values of all of the form Tui 

and therefore are equal to the corresponding tij variables. The theorem then follows from 
Corollary 14. 5[ allowing to substitute a min {n) for a max {n). □ 



4.3. A determinant formula for the generalized Lambda-determinant. To get a 

formula for the generalized Lambda-determinant, we simply have to combine the results 
of Theorems 14.41 and 14.61 Our generalized Lambda-determinant is the solution T 0i o,n of 
the T-system with coefficients subject to the initial conditions (12. 3p . It corresponds to the 
initial values: 

f 1 if z + j + n = mod 2 

otherwise 

V 2 > 2 

We assume from now on that the initial data t{j is chosen as in (14.31) . With the same 
notations as in the previous section, we obtain: 

Theorem 4.7. The generalized Lambda- determinant of any matrix A = {ai,j)i<i,j<n reads: 

n n 

\A\\ )I1 = a min (n) Y\ <\i Y\_ 

i=2 i=l 

As such, it is a Laurent polynomial of the aij's with coefficients in Z + [(Aj)jej, (fH)iej], 
where J — {1 — n, 2 — n, n — 1}. 

Proof. By direct application of Theorems 14.41 and 14.61 with the initial conditions (14. 3ft . As 
a>min( n ) is a principal minor of a product of Ui,U[,Vi,Vl matrices, and as the entries of 
all the matrices are Laurent monomials of the tj,/s and monomials of the Aj's and /i/s, 
we recover the Laurent polynomiality property of cluster algebras with coefficients, namely 
that the solution To,o,n is a Laurent polynomial of the entries of A, with coefficients that are 
polynomials of the Aj's and /Xj's. Moreover, from the form of the matrices, these polynomial 
coefficients have themselves non-negative integers coefficients. The theorem follows. □ 



5. From networks to 6V model 

5.1. Network formulation. The matrices V, V, U, U' may be interpreted as elementary 
"chips" that may be used to build networks, i.e. directed graphs with weighted edges. Each 
non-zero matrix element m it j of V, V', U, U' is interpreted as a the weight of a directed edge 
connecting an entry vertex i to an exit vertex j. This gives rise to the four following 
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"chips": 




where we have represented in dashed black lines the edges with weight 1, and in thickened 
red (resp. blue) edges the diagonals (resp. horizontals) carrying an extra multiplicative 
weight A (resp. fi). The SL r embedding Mj is clear: we place the chip for M in positions 
i, i + 1 and complete the graph by horizontal edges j — > j with weight 1, for j ^ i, The 
product of two matrices correspond in this language to the concatenation of the graphs, 
namely the identification of the exit points of the left matrix graph with the entry points of 
the right matrix graph. Graphs obtained by concatenation of elementary chips are called 
networks. Note that in the present case networks have face labels, that determine all the 
edge weights via the rules (15. ip . 

The matrix of the network is the product of the matrices of the chips forming it. The 
matrix element (i, j) of this product is the sum over all directed paths on the network with 
entry i and exit j of the product of the weights of the traversed edges. 

We may now represent the networks G m i n (n) and G max {n) corresponding to the matrices 
@mm(^) and m(1I (n). For n = 5, we have (recall that dashed black lines have weight 1, 
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while thickened blue (resp. red) lines have an extra weight \i (resp. A)): 
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G m in (5) 




and 




We have: 

Lemma 5.1. The function a m i n (n) (resp. a max (n) ) is the partition function for families of 
non-intersecting paths that start at the entry points 1, 2, n — 1 and end at the exit points 
1,2, ra-l of the network G min (n) (resp. G max {n)). 

Proof. By direct application of the Lindstrom-Gessel-Viennot Theorem [HI [12], as a min (n) 
(resp. a max {n)) are the n — 1 x n — 1 principal minors of Q m in(n) (resp. Q ma x(n)), the 
matrices of the networks G min {n) (resp. G max {n)). □ 

Remark 5.2. Note that due to the structure of G max (n) , there is a unique configuration 
of n — 1 non-intersecting paths contributing to a max (n), namely that in which all the paths 
have only horizontal steps. This gives an alternative pictorial interpretation of the proof of 
Theorem \4-6] 
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In view of Theorem 14.71 let us now reformulate a m i n (n). To this effect, we will simplify 
slightly the network G min {n) by use of the following: 

Lemma 5.3. We have the following equivalence between networks corresponding to the 
matrix product C/j_i(o, b, c; 1, l)Vi(b, c, d; A, fi): 

Xa/d 

i+ l o----jd - - -/>•• - - - ^ ac/bd 





where we have indicated the new edge weights (and kept the color code for thickened edges 
receiving extra multiplicative weights \x (blue) and A (red). 



Applying Lemma [5731 to all the pieces corresponding to products Ui-xVi in the network 
G m i n (n) allows to transform it into a new network G m i n (n) with the same network matrix, 
but which is now a subset of the directed triangular lattice (with edges oriented from left 
to right). For n = 4 this gives: 



1,4 



G m in (4) 




101 



where we have indicated the face variables corresponding to the initial data (I2.3I) . Note 
that the face labels of all up-pointing triangles are 1, whereas the down-pointing ones 
are the matrix elements of A. Note also that for technical reasons we have included a 
bottom isolated vertex. The matrix associated with this network, Q m i n (n), has therefore 
size 2n — 1 x 2n — 1. The entry point 1 is trivially identified with the exit point 1. 
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The edge weights read as follows: 



(5.2) 





(a) 



(b) 



depending on whether the triangle with face label a^- is on the SW border (a) for j — 1 
or anywhere else (b) for j > 1. In view of the formula of Theorem 14. 7\ it is natural to 
absorb the prefactor YYi=2 a Ti a redefinition of the network G min . Indeed, dividing the 
weights of the thickened colored edges in (a) by a i+ i t i, for i — 1, 2, n — 1 exactly absorbs 
the prefactor, as the paths starting at the left vertex must either traverse the red or the 
blue adjacent edge. After this transformation, all the weights have the form (b) of (15. 2p for 
1 < i, j < n — 1, including the ones on the SW border. Denoting by G m i n (n) the network 
thus modified, and by a m i n (n) the n x n principal minor of the associated matrix Q m i n (n), 
we have the following reformulation of the generalized Lambda-determinant: 

Theorem 5.4. The generalized Lambda- determinant of A is given by: 



\A\ 



X41 



n 



i=l 



As such, it is the partition function of families of n non-intersecting paths on the (directed 
triangular lattice) network G min {n), starting from the n SW vertices and ending at the n 
SE vertices, multiplied by all the entries of the last row of A. 

5.2. Another determinant formula for the generalized Lambda-determinant. 

Let us now briefly describe a method for explicitly computing the generalized Lambda- 
determinant of a given matrix A. The idea is to generate the paths eventually contributing 



to the minor a r 



n 



We consider points (x, y) in the integer plane Z , and a Hilbert space 



H with a distinguished canonical basis \x, y), x, y e Z, and dual basis (u, v\, u, v 6 Z such 
that (u, v\x, y) = 5 UiX 5 Vty . Let us introduce operators z,w acting on H as follows: 



(5.3) 



w\x, y) — \x + 1, y + 1) z\x, y) = \x + 1, y — 1) 



Consider now the embedding of the network G m i n (n) into Z 2 , with vertices (x, y) G D n , 
with: 



D, 



{(x,y) e Z 2 suchthatx + y = 0mod2and |x| + \y — n + 1| < n — 1} 
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Note that the bottom vertex of G m i n (n) is now at the origin (0, 0). The starting points of 
the paths contributing to a m i n (n) are Sj = (1 — i, i — 1), i — l,2,...,n and the endpoints 
are tj = (j — 1, j — 1), j = 1,2, ...,n. To each vertex (x, y) G D n , we attach the label of 
the adjacent down-pointing triangle of G m i n (n), {(x — 1, y + 1), (x + 1, y + 1), (x, y)}. For 
simplicity, we denote by a XjV = ax-y +n x+y +1 this label. We also attach the parameters 
£ y = \y + i- n and m x = n x respectively corresponding to the weights of the adjacent edges 
(x, y) — (x + 1, y + 1) and (x — 1, y + 1) — (x + 1, y + 1) of G m i n {n). This leads to the natural 
definition of the following diagonal operators: 

6c\x,y) = a Xj y\x,y), £\x, y) = £ y \x, y), m\x } y) = m x \x } y) 

Starting from a vertex (x, y) there are three possibilities of steps in G m i n (n), with the 
weights described in (15.21) (b). These weighted steps are generated by the following operator 
acting on H: 

Definition 5.5. We define the transfer operator T , acting on H as follows: 
T = w£ + wma~ 1 za + a~ 1 za 

T\x,y) = £ y \x + l,y + l) + m x+l \x + 2,y)+ \x + l,y-l) 

ot x +i,y-i a x+ i tV -i 

Let Zij(n) denote the partition function for paths from Sj to Cj on G m i n (n). 
Lemma 5.6. We have: 

Zij(n) = (j - T) _1 |l -i,i-l) 

or equivalently: 

Z id (n) = (0, Olw 1 -^/ - T)- 1 ^^, 0) 

Proof. By matching the Definition 15.51 with the weights of (15.21) (b), we easily see that each 
step of a path on G min {n) is generated by T. More precisely, the quantity (u, v \T k \x, y) 
is the sum over all paths on G min (n) of length k starting at the point (x,y) and ending 
at the point (u,v), of the product of the corresponding network edge weights. The first 
part of the Lemma follows. For the second part, we note that the adjoint = w" 1 , and 
w\x, y) = \x + 1, y + 1), so the dual is 

(x + 1, y + 1| = (x, y\w^ = (x, y\w~ l 

The second formula follows. □ 

Finally, combining Lemma 15.61 and Theorem I5.4[ we obtain a compact expression for the 
generalized Lambda-determinant: 

Theorem 5.7. 

(5.4) |^4|au= det {<d^\w l ~ j a(l -w£-wma~ 1 za-a~ 1 z&) z 1_i |0,0) 
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Proof. By definition, we have 

n 

a m in{n) = det Z itj (n) = )^ sgn(a) TT^rO;)j( n ) 

The formula of Theorem 15.41 includes a prefactor rij=i a ™,r Let us absorb the terms 
Oj_i,j-i = a nj - in each of the factors Z^^n), j = 1, 2, n. As (0, Olw) 1 "- 7 = (j — 1, j — 1|, 
we have <x,_ij_i(0, Olw 1 ^ = (0, Olw^a. The theorem follows. □ 

Remark 5.8. As a non-trivial check of Theorem \5. % let us evaluate the determinant in 
the homogeneous case when X a = X and \i a = fi for all a, and a^j = 1 for all i,j. All 
the operators z,w,£,rh commute, and b — I. Z it j(n) is therefore simply the coefficient of 
z w 3 of the series: 

i 

fz{z,w) 



1 — z — Xw — fizw 

where the subscript Z stands for the (infinite) matrix with entries Zjj(oo). Note that fz is 
independent of n, which is simply the size of the principal minor of Z we are supposed to 
compute to obtain \A\\^. We have the following (infinite size) LU decomposition: 

Z = B(1,1)*B(A + ^,A) 

where the infinite matrices B(a,B) are upper triangular, with non-zero matrix elements 
B(a,(3)i t j = Rja 4 /?- 7 ' - ' for j >i>0, or equivalently with generating functions 



fB(a,/3)(z,w) = ^ B ( a 'l3) 



;Z w 



M 1 — Bw — azw 

i,j>0 ' 

The LU decomposition is easily obtained by the convolution formula: 

/dt 

where the integral extracts the constant term in t. Finally, by triangularity of the factors, the 
product Z = B(l, I) 1 B(X + fi, X), when truncated to the first n rows and columns, truncates 
to the product of the n x n truncations of the factors. We deduce that the principal n x n 
minor is: 

14a* = det (5(1, 1U det (B(X + fx, X) itj ) = (A + /i)" ( "~ 1)/2 

0<t,3<n— 1 0<i,]<n— 1 



In agreement with the result of Example \3.4\for a% = 1 for all i. 
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5.3. Transformation into 6V model/ ASM: proof of Theorem 11.21 By Theorem 15.41 
the generalized Lambda-determinant is expressed as a sum over configurations of n non- 
intersecting paths on the network G m i n (n). Any such configuration determines entirely 
the state of the down-pointing triangles labelled by the matrix elements (ciij)ijei (we may 
draw such triangles around a^i, a^ n a%j and a n j as well). We say that a vertex (resp. 
edge) is occupied (resp. empty) if a path of the configuration traverses it. This gives rise 
to seven possible local configurations of the bottom vertex and the three edges of each 
down-pointing triangle, that we associate with the local configurations of the 6V model as 
follows: 




where we have represented occupied vertices/edges by filled circles/thick solid lines, and 
empty ones by empty circles/dashed lines. This is a mapping 2 m to 1 of the configurations of 
n non- intersecting paths on G m i n (n) with a total of m down-pointing triangle configurations 
among the two on the right to those of the 6V model with DWBC on the nxn grid rotated 
by +vr/4, with m vertices of type c<i- The latter in turn are in bijection with ASMs of size 
nxn with m entries equal to —1. This is illustrated in the case n = 4, m — 1 below: 




Let us now collect the various weights of the path configurations contributing to 
by distinguishing the dependence on the matrix elements of A (face labels) and the co- 
efficients A, n (edge labels). The matrix elements the labels of the down-pointing 
triangular faces. For each face labeled a, only four edges of the network G m i n (n) have 
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weights depending on a. These are: 




where we have shaded the down-pointing triangle with label a, and indicated the a- 
dependence of the four relevant edges. This gives rise to the only following contributions 
involving the matrix element a, depending on the state of occupation of the four above 
edges, which we list by configuration of the down-pointing triangle: 



(5.6) 



6V 




a 2 




b 2 




c 2 






\_ 






* 


\ a 

a 


network 


'/ 






\ a 

*"a**- 












\ a 

\ 




\ a 

\ 




weight 


1 


1 


1 


1 


a 


a" 1 



The dependence on the face label a = ay is therefore a b , where b = bij G {0, ±1} is the 
matrix element of the ASM in bijection with the associated 6V model configuration. This 
however is only valid for i ^ n, and extra care should be taken with the SE down-pointing 
triangles. Indeed, their bottom vertex must always be occupied (it is an exit point for the 
non- intersecting path configuration), and only two configurations may occur: (i) with the 
top edge occupied (6V configuration b\) or (ii) with the left diagonal edge occupied (6V 
configuration a 2 ). In both cases, the occupied edge receives a weight a -1 : 




However, recall that the formula for the generalized Lambda-determinant of Theorem 15.41 
contains a prefactor Y]a=i a n,i which is the product over the labels of the SE down-pointing 
triangles. We may therefore absorb this prefactor in a redefinition of the weights of hor- 
izontal and left diagonal edges adjacent to the SE faces, namely by multiplying them by 
the face label a. The resulting network G' min {n) and the n x n principal minor a' min {n) of 



2(> 



PHILIPPE DI FRANCESCO 



its matrix are such that 

(5.7) \A\ X ,» = a' min (n) . 

Note that this latter transformation of the network into G' min {n) corresponds exactly to the 
insertion of the a operator in Theorem 15. 71 According to the table (15.61) . each configuration 
of n non- intersecting paths on the network G' min {n) that contributes to a' min (n) receives a 
weight ■ a\y where bij is the ASM in bijection with the 6V configuration associated to 
the path configuration. 

Finally, the dependence of the edge weights on the coefficients \ a , /i a is clear from (15.21) 
(b): around the down-pointing triangle with label a^j, the horizontal edge receives the 
weight fjLi+j- n -i, while the left diagonal edge receives the weight Xj-i. This gives the total 
weights 



(5.8) 



6V 


»1 


a 2 


bi 


b 2 




c 2 


network 


V 


\ 


"\ a / 
V 


a' 


\ a 
*' 




weight 




1 


n 


1 


a 


^ia 1 A, a 1 



where for a = dij, X = Xj-i and fi = for all i,j = 1, 2, ...,n. 

Starting from the expression ( 15.71) . we may rearrange the configurations of paths with a 
total of m triangle configurations of type C2 say in positions (ii, ji), (im,jm) ( an d 

a total of m + n triangle configurations of type Ci in fixed positions) into a 2 m -tuple, in 
bijection with the unique ASM with exactly m entries —1 at the same positions (and m + n 
entries equal to 1 at the same positions). Summing over the weights for such a 2 m -tuple 
yields an overall factor of 

m 

\\ (^jr-ir + fMr+jr-n-l) 
r=l 

according to the table (15. 8p . Collecting all the other weights, we may rewrite \A\\ ;/jl as a 
sum over n x n ASMs B, with the weights of the form Yli j£i v i,j( a i,j^ij\^j-iy ^i+j-n-i) 
where 

r A if Q id (B) = a, 

, , i , s h I u if Cj AB) = b\ 
ViJ {a,b\\,n) = a x A + /i . f g ;2b) = c 2 

I 1 otherwise 

This completes the proof of Theorem 11.21 

6. Discussion and conclusion 

In this paper, we have defined the generalized Lambda-determinant of any square matrix 
A. Via Theorem ll.2l this quantity may be reexpressed in terms of weighted configurations 
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of the 6V-DWBC model on a square grid of the same size as A, in bijection with ASMs of 
the same size. 



6.1. Properties. Let us show how the formula of Theorem 11.21 may be applied to give an 
alternative proof of Proposition 13.51 Notations are as in Proposition 13.51 

We simply have to examine how the 6V/ASM configurations contributing to the formula 
(11. 2p are changed when we change A — > tp(A), namely (i) rotate the matrix A by a clockwise 
quarter-turn —tt/2 {ip = cr) or (ii) transpose the matrix A, namely reflect it w.r.t. to its 
diagonal (ip = r). We have the following action of ip on the 6V-DWBC configurations C, 
compatible with the transformation B — > ip(B) for any n x n ASM B: 

Lemma 6.1. Let if £ {a, r} act on configurations C of the n xn grid 6V-DWBC model 
as follows: <£>(C) obtained by applying directly the transformation <p to the arrow config- 
uration (rotation by —tt/2 for <p = a, reflection w.r.t. the diagonal for ip — r), and then 
performing a global flip of all the edge orientations. Then we have for all n x n ASM B: 



<p(e(B))=e(p(B)). 



Proof. Note first that the global sign flip reinstates the correct orientations of external 
edges to restore DWBC boundary conditions. Under the transformation <p, the six vertex 
configurations are mapped as follows: 

T (a 1 )=a 1 r(&i) = b 2 r(a 2 ) = a 2 r(b 2 ) = h r(c x ) = c x r(c 2 ) = c 2 

As both Ci and c 2 are invariant under (p £ {a, r}, the ASM associated to the configuration 
<p{G) is simply <p(B), and the Lemma follows. □ 



Let us now evaluate ^(^^(a),^)- Using the formula (II. 2p . the weight w reads: 

r vWj-i if e(B) id = ai 



)i+j—n—l 



ae(B) iJ = c 2 



1 otherwise 
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Let us perform the change of summation variables B = (p(C). For ip = o~,t respectively, 
this gives weights: 



These two expressions match Wij(A,C; A,/i) provided the 6V weights in the definition are 
changed into respectively (i) a\ — > a 2 and (ii) &i — > b 2 . 

Proposition 13.51 follows from the fact that the formula (11.21) remains valid if we change 
a\ — > a 2 and/or b\ — > b 2 in the definition of w. This in turn is the consequence of the 
following: 

Lemma 6.2. For any 6V-DWBC configuration G, we have for all a: 

\{i G {1, 2, n} such that C^+a = ai}\ = \{i G {1, 2, n) such that Qi t i +a = a 2 }\ 
\{i G {1, 2, n} such that Qi <a -j = ^i}| = |{* G {1, 2, n) such that Qi, a -j = &2j-| 

in other words there are as many a\ configurations as a 2 along each parallel to the diagonal, 
and there are as many b± configurations as b 2 along each parallel to the anti- diagonal. 

Proof. Let us consider a diagonal (i,i + a)i<j<„_ a of C, and the sub-configuration situated 
under the line (i, i + a — 1) (see Figfj] (a)). This configuration has entering edges on its left 
and possibly right vertical borders, and outgoing edges along its bottom and possibly top 
horizontal borders, in equal numbers. By the neutrality condition at each vertex (there are 
as many entering as outgoing edges) , we deduce that there are as many entering as outgoing 
edges along the diagonal border. These are the N and E edges of the vertices of the diagonal 





1 otherwise 




Equivalently, we may write: 
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J i L J i 1 1 1 1 1 1 1 L 

» » • • • • 1 — ». , a 

H !!!!!!. i_ t! ! 



• • • • 

• • • 



• • • r 

• • r 

• r 

r 

TTTTTTTTTr ^ I r"" 

(a) (b) 

FIGURE 1. The sub-configuration of a 6V-DWBC configuration C situated (a) 
under a diagonal (b) above an anti-diagonal. We have marked the edges of the 
diagonal border, respectively equal to the N and E edges (a) and S and E edges 
(b). 

border. Inspecting the vertices on ( 11. II) we see that all vertices have opposite orientations 
of their N and E edges, except for ai and a 2 who have respectively two outgoing and two 
entering N and E edges. We deduce that these occur in the same number to guarantee 
that as many edges are entering and outgoing along the diagonal border. The argument 
may be repeated for anti-diagonals by considering the sub-configuration situated above the 
anti-diagonal (i, a + 1 — i) (see Fig{T] (b)). In that case, we must inspect the S and E edges 
of the vertices of (II. ip . and we note that only b\ and b 2 have S and E edges pointing both 
out or in, therefore they must be in equal number along the diagonal border. The Lemma 
follows. □ 

6.2. The generalized Lambda-determinant in terms of ASM. Theorem 11.21 ex- 
presses the generalized Lambda-determinant in terms of 6V-DWBC configurations. Let 
us reexpress this quantity purely in ASM terms. To this end, we must better explain 
the inversion number of any ASM B. The contributions to Inv(5) = X^<j k<£bj,kh,e 

\k C\ 

are obtained from rectangular submatrices of B^'J of B, obtained by retaining only the 
intersection between rows + 1, j) and columns (k, k + 1, ...,£) in B. To produce a 
non-trivial contribution to Inv(5) the rectangle must have non-zero elements in position 
(j, k) (SW) and (i, €) (NE). Let us now concentrate on the element in position (z, k) (NW). 
Let h, v denote respectively the first non-zero element of B on row i to the right of 6^ and 



30 



PHILIPPE DI FRANCESCO 



on column k below 6^. The simplest way of thinking about the bijection from ASM to 
6V-DWBC is that along rows edges are oriented from the entries —1 to the entries 1 and get 
reversed at the first and last 1, with opposite rules along columns. Assuming = 0, this 
gives the following 4 possibilities: (i) (h, v) = (1, 1): then Q(B) i k = d; (ii) (h,v) = (—1, 1): 
then G(B) itk = b 2 ; (iii) (h,v) = (1,-1): then 6(5)^ = h; (iv) (h,v) = (-1,-1): then 
G(B)i t k = ci2- When computing contributions to Inv(S), we may first fix the NW corner 
of rectangles we sum over, and finally sum over the position (z, k) of that corner. In cases 
(ii-iv), for fixed (i, k) the sum gives trivially zero, as there are as many l's and —1 either 
on the row or the column of (i, k) or both. So only the case (i) contributes a total of +1 
to Inv(.B). Moreover, if b^ ^ 0, only cases (i) and (iv) survive, with respectively b it k = — 1 
and bi t k = 1. However for the same reason as before, the case (iv) does not contribute to 
lrw(B) when we sum over the rectangles with fixed NW corner b^k- We are left with only 
the case (i), for which Q(B)ij = c 2 , and the sum over rectangles contributes a total of +1 
to inv(B). 

To summarize, terms contributing +1 to Inv(B) may be associated bijectively to the a\ 
and c 2 configurations in C(B). So we have lnv(B) — #(— l)s = #(oi), the total number 
of ai vertices in G(B). Moreover, the vertices equal to a± in G(B) correspond bijectively to 
the entries of B, such that the first non-zero elements in B in their row to the right and 
in their column below are both 1. 

We may repeat the above analysis with the matrix cr(B), rotated by — n/2. What plays 
the role of ai after rotation are the vertices b± before rotation (see previous subsection), and 
C2's remain unchanged. We deduce that lnv(a(B)) — #(— 1)b = #(&i), the total number 
of b\ vertices in G(B). Moreover the vertices equal to b± in Q(B) correspond bijectively to 
the entries of B, such that the first non-zero elements in B in their row to the right and 
in their column above are both 1. 

Definition 6.3. For each fixed diagonal 5 a (B) = {b it j | j — i = a} for 2 — n < a < n — 2, 

let us denote by l a {B) the total number of entries such that the first non-zero entries 
to the right and below are both 1. For each fixed anti-diagonal otb = {hj \ i + j '• = b} for 
2 — n<b<n — 2let l' b (B) denote the total number of entries such that the first non-zero 
entries to the right and above are both 1. Finally, let #(— 1)5(2', j) = 

We have the following generalized Lambda-determinant expression: 

Theorem 6.4. The generalized Lambda- determinant of any n x n matrix A reads: 

= e n xi a [B) ti iB) n * (v* + ^-,-i) #( - i)s(ij) 

nxnASMB a=2-n M=l 

6.3. The general solution of the T-system with coefficients. In this paper, we have 
concentrated on special initial conditions for the T-system (13.21) . namely such that Tj ^o = 1 
and Ti Starting from the general T-system solution of Theorem I4.6[ 
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where both and T^i are assigned arbitrary initial values Uj = Ti t j jCi j , we may 
repeat the steps we took to modify the initial network G m i n {n) into a subset of the directed 
triangular lattice. The resulting network (which we still denote by G m i n {n), with associated 
minor a min {n) by a slight abuse of notation) has the following edge weights around each 
up-pointing triangle: 



(6.2) 



and we have 




where the last transformation to G' min {n) consists in absorbing the tj- n j-i in the SE hori- 
zontal and diagonal edge weights adjacent to the endpoints of the paths. Noting that the 
initial values 7^,1 are the face labels of the down-pointing triangles, while the Tij t o are 
those of the up-pointing triangles, we may summarize the dependence of the edge weights 
on the face labels of the down-pointing and up-pointing triangles as follows: 




where a,b — tij for some Combining these the weights of down-pointing triangles (15. 6p . 
we note that the weights acquire a simpler form if we transfer the A and /i parameters from 
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the weights of down-pointing triangles to those of up-pointing ones, resulting in: 



network 



(6.3) 



weight 1 
network 
weight 1 



V 



a 



a 



a 



a 

V 



a 



a 

V 



a 



^a 



a 



A, a 



-1 



where for a 



£jj, we have A 



Xi and /i = fij. 



Remark 6.5. We have associated a 6V-DWBC configuration of size n to each configuration 
of n non-intersecting paths on G' min {n) from SW to SE, by the local rules ( 15.5ft associating 
6V vertices to each down-pointing triangle configuration. We could have done the same for 
up-pointing triangles as well, by using the same rules on the reflected up-pointing triangles 
w.r.t. their horizontal edge. More precisely, this would yield a 6V-DWBC configuration 
of size n — 1 where the prime refers to the fact that the boundary condition has opposite 
convention (horizontal external edges point out of the grid, while vertical ones point into 
the grid). Flipping all edge orientations, this becomes a 6V-DWBC configuration of size 
n — 1. So each configuration of non-intersecting paths on G' min (n) that contribute to a ; min (n) 
give rise naturally to a pair (B,C) of ASMs of respective sizes n and n — 1. Such pairs 
are called compatible pairs in the literature. However, these pairs do not play any role in 
our fully inhomogeneous T -system solution, as weights are attached only to such compatible 
pairs and cannot be disentangled nicely. In particular there is no resummation of the non- 
intersecting lattice path weights into weights only pertaining to the ASMs, as opposed to 
the case where T^o = 1 of Theorem ! where the second ASM plays no role. 



Example 6.6. We consider n 
the plane): 



3. We have the initial data assignments (represented in 



To,- 



-2,1 



t 



0,-2 



^1,-1,1 

To,-i,o 



^2,0,1 

£1,-1 ^1,0,0 

^0,-1 ^0,0,1 
t-1,-1 ^-1,0,0 
^-2,0,1 



^2,0 
tl,0 

to,o 

t-1,0 
t-2,0 



^0,1,0 



^0,1 



To, 



2,1 



1-0,2 



There are 8 distinct configurations of the 3 non-intersecting paths on G' min (3), contributing 
to a' min (3) (we have indicated the by small squares the centers of the n 2 = 9 down-pointing 
triangles and of the (n — l) 2 = 4 up-pointing ones, corresponding to labels tij with (i,j) G 
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=4 



-2-10 1 2 

(a) 



^•3 



-• 2 



(b) 



p-- 

3»- -4 * 

2* ^2 



(c) 



2» *2 

1» 



(d) 



3 -p;-' ^^Ni ; 





(e) 



(f) 



(g) 



(h) 



T/ie corresponding weights read: 



„,, — ,. ,. ,. *0,-2.to,0,*0,2 

W(a) - V-lVoVl t0 _ lt01 

„„ — \2,. -l*-l,l*l, -1*1,1 

W W - A 0/ i i ,-i*o,i*o,o 

*l,-l*l,l*-2,0 



4-1,0*1,0 



W(6) = A /i /i; 
lt>( e ) = A /io 



*-!,-!*!, -1*0, 2 

*0, -1*0,1 
2 t-l,-lt-i, 1*1, -1*1,1 



*-l, 0*1, 0*0,0 
*-!,-!*-!, l*-2,0 

*-l, 0*1,0 



W(c) = /U-l/UoAo 

W(j) = a* AoAi 



*0,-l*-l, 1*1,1 

*0, -1*0,1 
*-l,-l*-l, 1*2,0 

t-1, 0*1,0 



The corresponding solution of the T -system is: 



rp to,-2, ^0,0; t-0.2 , ^-1,-1^1,-1^0,2 x ^0,-1^-1,1^1,1 

J- 0,0,3 — p*— lMo/J-l T 7 r Ao/^o/^l 7 7 r /X-i/^o^O 



+A /io 



^0,-1*0,1 

£-i, —l^—i, A, 

^0,-1^0,1^0,0 



+ AoyUo 



^0,-1^0,1 
2 £-1, -1^-1,1^1, -1^1,1 



-A_iAoA*o — L ; ^ : — I" A-iAqAi 



£-1,0^1,0^0,0 
£-1,-1^-1, 1^-2,0 



AtoAoA 



to, -1^0,1 

£-i,-i£-i,i£2,o 
£-i,o£i,o 



1,0 



£-1,0^1,0 
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If we take t±i t o = £o,±i = 1 ; we recover the formula of Theorem \1.2t 



,0,3 — A*— lA*0A*l ^0,-2) ^0,0? ^0,2 + ^OpOpl ^-1,-1^1,-1^0,2 + /^-l/^O-^O ^0,-1^-1,1^1,1 

+/ioAoAi t-\ t —1^—1,1^2,0 + A_iAo)Uo ^i, -1^1,1^-2,0 + A_iAqAi t— 1,— it— 1, it-2,0 




where the sum extends over the 7 ASM 0/ sue 3. Finally, for tij = 1 for all i,j and 



in agreement with Theorem \3.6\ for n = 3. 

6.4. Conclusion and perspectives. A number of questions remain to be answered con- 
cerning this new deformation of the determinant. 

The partition function of the 6V-DWBC model at its free fermion point A = 0, q = e"^ 2 , 
corresponds to the 2-enumeration of ASM, with a weight 2 per entry —1, and matches the 
Lambda-determinant, with all A a = Pb — 1. However, the 6V-DWBC model is known to 
have a non-trivial deformation involving 2n spectral parameters say Zi per row i and Wj per 
column j of the square lattice nxn grid. Our generalized Lambda-determinant introduces 
2n — 2 other deformation parameters, which are attached to diagonals and anti-diagonals 
of the same square grid. It would be interesting to mix the two deformations and see if we 
can obtain some new information on ASMs by the process. 

Another question concerns the limit shape. It is known that the 6V-DWBC model at 
its free fermion point is also in bijection with the domino tilings of the Aztec diamond [8], 
for which an arctic circle theorem exists, namely there exists a limiting curve (a circle) 
separating frozen tiling phases from entropic tiling in the limit of large grid and small mesh 
size. It is easy to work out the homogeneous deformation of this curve when A a = A and 
Pb — p for all a, b. The simplest way of computing the artic curve in our model (in the 
case of a matrix A with entries a^j = 1) is by analyzing the singularities of the generating 
function of the density pij^ defined as 



where Tij^ is the solution of the T-system (13 .2D with initial data (14. 2ft . Defining p(X, Y, Z) = 
J2k>i i j&Pi,j,kX % Y 3 Z k we find, by differentiating the T-system relation and evaluating the 



Aj = fii — q l we have: 



T 0)0 ,3 = 4 + 2(q + q- 1 ) = (1 + ^(l + q°)(l + q 1 ) 



Pi,j,k 



dLogT itjtk 



result at L 



h3 ~ 



1: 



p(X,Y,Z) 
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Following [13], we find that in the coordinates (x, y) obtained as limit of (|, |), for (i, j) e 
Dij jk , of a kxk 6V-DWBC grid rotated by +vr/4, the singular curve is the following ellipse: 



x \ 1 + -,) + v'{ 1 + J 

It would be interesting to see how this is affected by inhomogeneous parameters A a ,/i a . 
In the particular example \ — pi — q l of Section 13.31 we find the following equation for 
p(X,Y,Z): 

(1 + Z 2 -X- X- 1 )p(X, Y, Z) + (l + Z 2 -qY- q' x Y '~ x ) pig' 1 X , qY, Z) = 2Z 

For generic q (not a root of unity) the function p(X, Y, Z) is not algebraic, and may be 
defined as the following limit: 

n(X YZ)- 22 lim fV-rt* 17 O^ygKWZZ^g) 

P[ ' ' J " (l-XZ)(l-X-!Z)t™-^ 1 J ll(i_ ? -m 

It would be interesting to investigate its singularity structure. 

Finally, the connection to cluster algebra should also allow us to consider quantum 
versions of the T-system, via the quantum cluster algebra construction of [T]. A first step 
in this direction was taken in [5] for the case of the Ai quantum T-system. Presumably, 
such quantum deformations should yield new non-trivial deformations of the Lambda- 
determinant and hopefully enrich our understanding of ASMs. 
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